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ABSTRACT 

In mammalian mitochondria, 22 species of tRNAs en- 
coded in mitochondrial DNA play crucial roles in the 
translation of 13 essential subunits of the respiratory 
chain complexes involved in oxidative phosphoryla- 
tion. Following transcription, mitochondrial tRNAs 
are modified by nuclear-encoded tRNA-modifying 
enzymes. These modifications are required for the 
proper functioning of mitochondrial tRNAs (mt tR- 
NAs), and the absence of these modifications can 
cause pathological consequences. To date, how- 
ever, the information available about these modifi- 
cations has been incomplete. To address this is- 
sue, we isolated all 22 species of mt tRNAs from 
bovine liver and comprehensively determined the 
post-transcriptional modifications in each tRNA by 
mass spectrometry. Here, we describe the primary 
structures with post-transcriptional modifications of 
seven species of mt tRNAs which were previously 
uncharacterized, and provide revised information re- 
garding base modifications in five other mt tRNAs. 
In the complete set of bovine mt tRNAs, we found 
15 species of modified nucleosides at 118 positions 
(7.48% of total bases). This result provides insight 
into the molecular mechanisms underlying the de- 
coding system in mammalian mitochondria and en- 
ables prediction of candidate tRNA-modifying en- 
zymes responsible for each modification of mt tR- 
NAs. 

INTRODUCTION 

Mitochondria, organelles present in most eukaryotic cells, 
provide the chemical energy required by living cells in the 
form of adenosine triphosphate (ATP), which is synthesized 
by the electron transport chain and oxidative phosphoryla- 
tion (1). Mitochondria contain their own genomic DNA, 
called mitochondrial (mt)DNA, and unique transcription 
and translation machinery that converts their genetic infor- 



mation into proteins. In mammals, the mtDNA is a circular 
double-stranded DNA, ~16 kilobase pairs (kb) in length, 
which contains 1 3 genes encoding essential subunits of the 
respiratory chain complexes and 24 RNA genes (2 riboso- 
mal RNAs and 22 tRNAs) required for mitochondrial pro- 
tein synthesis. 

The mammalian mitochondrial decoding system differs 
from the canonical decoding system by its use of four non- 
universal codons (Table 1) (2): AUA for Met, UGA for Trp 
and AGR (R = A or G) for Stop. The 60-sense codons are 
deciphered by 22 species of mitochondrial tRNAs, which 
constitute the smallest set of tRNAs necessary to trans- 
late all sense codons among all kingdoms of life, includ- 
ing other organelle decoding systems. Post-transcriptional 
modifications at the first letters of tRNA anticodons play 
a critical role in establishing this minimal decoding sys- 
tem. To reduce the number of tRNA species, each of eight 
family boxes in mitochondria is decoded by only a single 
tRNA. The tRNAs responsible for the family boxes fre- 
quently have unmodified uridines (U34) at the first (wob- 
ble) position of an anticodon. According to Crick's wob- 
ble rule (3), U34 can recognize only A and G at the third 
position of a codon. In the decoding systems of some bac- 
teria and most organelles, U34 can base-pair with any of 
the four bases by the so-called 'four-way wobble rule' (2) 
or 'super wobbling' (4). The conformational flexibility of 
U34 is thought to enable U:U and U:C pairing. On the 
other hand, tRNAs responsible for two-codon sets ending 
in purines (NNR; N = any four nucleotides) have modi- 
fied uridines (xm 5 s 2 U-type) at their wobble positions (5). 
In general, xm 5 s 2 U-type modification restricts conforma- 
tional flexibility of the wobble base, thereby strengthening 
recognition of NNR codons and preventing misrecognition 
of codons ending in pyrimidines (NNY; Y = U and C). 
We previously discovered that taurine-containing modified 
uridines in mammalian mt tRNAs are responsible for NNR 
codon sets: 5-taurinomethyluridine (Tm 5 U) at the wobble 
position of human and bovine mt tRNA Leu ^ UUR ^ (6), and 
5-taurinomethyl-2-thiouridine (Tm 5 s 2 U) at the wobble po- 
sition of human and bovine mt tRNA Lys (6). Subsequently, 
Tm 5 s 2 U was also found at the wobble position of the bovine 
mt tRNAs for Glu and Gin (7). Determining the proper us- 
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age of unmodified or modified U34 in tRNAs would reveal 
a fundamental principle of the minimal decoding system in 
mammalian mitochondria, nullnull 

5-formylcytidine (f^C), another unique modification in 
mammalian mitochondria (8), is present at the wobble po- 
sition of mt tRNA Met . Biochemical studies show that f^C is 
required for recognition of the non-universal AUA codon, 
in addition to the canonical AUG codon (9,10). Possi- 
ble base pairing between f^C and A was demonstrated 
by a crystallographic study using an f^C-containing anti- 
codon stem-loop bound to the 30S subunit of Thermus ther- 
mophilus (11). 

Defective mitochondrial translation results in mito- 
chondrial dysfunction, ultimately causing pathological 
consequences (12). Numerous pathogenic mutations as- 
sociated with mitochondrial diseases have been found 
in mtDNA (http://www.mitomap.org/MITOMAP) (13). 
These pathogenic mutations are maternally inherited. 
Among over 400 pathogenic mutations compiled to date, 
~200 have been mapped to mt tRNA genes; these mu- 
tations diminish biogenesis, stability and function of tR- 
NAs (2,12,14-15). Two major subgroups of mitochondrial 
encephalomyopathies are caused primarily by point mu- 
tations in mt RNAs: mitochondrial encephalopathy, lac- 
tic acidosis and stroke-like syndrome (MELAS), caused by 
a mutation in mt tRNA Leu ( UUR ); and myoclonus epilepsy 
with ragged-red fibers (MERRF), caused by a mutation 
in mt tRNA Lys . Approximately 80% of MELAS patients 
have an A-to-G mutation at position 3243 (A3243G) (16) 
in the mt tRNA Leu(UUR) gene, and another 10% have a 
T-to-C mutation at position 3271 (17). MERRF patients 
have an A8344G mutation in the mt tRNA Lys gene (18). 
We previously reported that im 5 U and Tm 5 s 2 U were not 
present in mutant mt tRNA Leu(UUR) harboring the A3243G 
or T3271C mutation (19), or in mt tRNA Lys harboring the 
A8344G mutation (20). We also confirmed the absence of 
taurine modifications in MELAS patient tissues harboring 
A3243G, G3244A, T3258C, T3271C or T3291C mutations 
(21), as well as in MERRF patients harboring the A8344G 
mutation (22). These pathogenic point mutations are as- 
sumed to act as negative determinants of im 5 (s 2 )U biogen- 
esis. In MELAS, the absence of im 5 U in mt t RNA Leu(UUR) 
results in a defect in decoding the UUG codon, leading to 
lower expression of the UUG-rich protein ND6 (23). Simi- 
larly, in MERRF the absence of both 5-taurinomethyl and 
2-thio groups of Tm 5 s 2 U in mt tRNA Lys leads to severe 
translation failure of both types of AAR codon (24). These 
observations imply that deficiency in modification of mt 
tRNA plays a key role in molecular pathogenesis. As men- 
tioned earlier, unmodified U34 can read any of four bases 
at the third letter of codons in a family box by the four-way 
wobbling. However, unmodified U34 is only used in tRNAs 
responsible for family box codons in which at least one G 
or C is present at the first or second letter of codons. If the 
codon-anticodon interaction is stabilized by one or two GC 
pairing at the first two base pairings, U34 is considered to 
read any of four bases in the family boxes. In the case of 
mt tRNA Leu(UUR) and mt tRNA Lys , there is no G or C at 
the first or second letter of their cognate codons. This is the 
reason why unmodified U34 in the mutant tRNAs cannot 



decipher cognate codons efficiently, not by expanding their 
decoding capacity. 

Mitochondrial diseases are also caused by pathogenic 
mutations in nuclear-encoded genes (2), including genes en- 
coding translation factors, aminoacyl-tRNA synthetases, 
tRNA processing enzymes and tRNA-modifying enzymes. 
Loss-of-function mutations in these genes hamper the bio- 
genesis and function of mt tRNAs. Several instances of 
pathogenic point mutants in tRNA-modifying enzymes 
have been reported to date: mitochondrial myopathy and 
sideroblastic anemia (MLASA), acute infantile liver failure 
and hypertrophic cardiomyopathy with lactic acidosis are 
associated with pathogenic mutations in PUS1 (25), MTU1 
(26) and MTOl (27,28), respectively. 

To gain more insight into the molecular basis of the 
mitochondrial decoding system, as well as the molecular 
pathogenesis of human diseases caused by deficiencies in 
mt tRNA modifications, it is necessary to obtain a com- 
plete picture of post-transcriptional modifications in all 22 
species of mammalian mt tRNAs. To date, 1 1 species of hu- 
man or bovine mt tRNAs have been sequenced and their 
post-transcriptional modifications determined (2). In pre- 
vious work, we isolated all 22 species of bovine mt tRNAs 
(29). By analyzing these tRNAs by mass spectrometry, we 
determined the post-transcriptional modifications of seven 
species of mitochondrial tRNAs that have never been re- 
ported, and provide some revised information regarding the 
modified bases in five other mitochondrial tRNAs. In to- 
tal, we identified 15 species of modified nucleosides at 118 
positions in 22 species of bovine mt tRNAs. We discuss 
the basic principles of the mitochondrial decoding system 
in mammals, and propose candidate tRNA-modifying en- 
zymes whose roles remain to be confirmed experimentally. 

MATERIALS AND METHODS 

Isolation of individual mitochondrial tRNAs from bovine liver 

Bovine liver RNA was prepared as described previously 
(6,29). Briefly, crude RNA was extracted from buffer- 
homogenized bovine liver by phenol extraction, and the 
tRNA fraction was roughly concentrated by anion ex- 
change chromatography. Individual tRNAs were isolated 
by chaplet column chromatography (29). Seventeen of them 
were isolated homogeneously based on the polyacrylamide 
gel electrophoresis analysis. Five mt tRNAs for Leu(UUR), 
Asn, Thr, Met and Pro were isolated as major bands 
but had some minor bands. However, the qualities of 
the isolated tRNAs are sufficient enough to analyze post- 
transcriptional modifications. 

Cyanoethylation of pseudouridine in tRNA 

Cyanoethylation of tRNA was performed basically as de- 
scribed (30). Eight micrograms of isolated tRNA dissolved 
in less than 4 (xl of Milli-Q water was added to 30 |xl of 41% 
(v/v) ethanol/1.1 M trimethylammonium acetate (pH 8.6). 
After addition of 4 (xl of acrylonitrile (Wako Pure Chemi- 
cal Industries), the mixture was incubated at 70°C for 2 h, 
lyophilized and dissolved in Milli-Q water. The solution was 
then subjected to RNase digestion and analyzed by RNA 
mass spectrometry, as described below. 



7348 Nucleic Acids Research, 2014, Vol 42, No. 11 



Table 1. Codon-anticodon pairing in the bovine mitochondrial genetic code 



Codon 


Amino acid 
(anticodon) 


Codon 


Amino acid 
(anticodon) 


Codon 


Amino acid 
(anticodon) 


Codon 


Amino acid 
(anticodon) 


UUU 


Phe 


ucu 




UAU 


Tyr 


UGU 


Cys 


uuc 


(GAA) 


ucc 


Ser 


UAC 


(QUA) 


UGC 


(GCA) 


UUA 


Leu 


UCA 


(UGA) 


UAA 


stop 


UGA 


Trp 


UUG 


(WUAA) 


UCG 




UAG 




UGG 


(Tm 5 UCA) 


CUU 




ecu 




CAU 


His 


CGU 




cue 


Leu 


ccc 


Pro 


CAC 


(QUG) 


CGC 


Arg 


CUA 


(UAG) 


CCA 


(UGG) 


CAA 


Gin 


CGA 


(UCG) 


CUG 




CCG 




CAG 


(tiuVUUG) 


CGG 




AUU 


He 


ACU 




AAU 


Asn 


AGU 


Ser 


AUC 


(GAU) 


ACC 


Thr 


AAC 


(QUU) 


AGC 


(GCU) 


AUA 


Met 


ACA 


(UGU) 


AAA 


Lys 


AGA 




AUG 


(f 5 CAU) 


ACG 




AAG 


(tiuVUUU) 


AGG 


stop 


GUU 




GCU 




GAU 


Asp 


GGU 




GUC 


Val 


GCC 


Ala 


GAC 


(QUC) 


GGC 


Gly 


GUA 


(UAC) 


GCA 


(UGC) 


GAA 


Glu 


GGA 


(UCC) 


GUG 




GCG 




GAG 


(TmVUUC) 


GGG 





Non-universal genetic codes are denoted in bold type. AUA: He (universal), Met (mitochondria); UGA: stop (universal), Trp (mitochondria); AGA/G: 
Arg (universal), stop (mitochondria). 



RNA mass spectrometry 

RNA fragments digested by RNases were analyzed by mass 
spectrometry as described previously (31,32). In brief, 2-5 
ng of isolated tRNA was digested with RNase Ti (Epicen- 
tre) or RNase A (Ambion) and analyzed by an LTQ Or- 
bitrap mass spectrometer (Thermo Scientific) with a nano- 
electrosprayer connected with a splitless nanoflow high 
pressure liquid chromatography system (DiNa, KYA Tech- 
nologies). Alternatively, 0.4-2 (xg isolated tRNA was di- 
gested with RNase Ti and analyzed on a LCQ DUO ion- 
trap mass spectrometer with an ESI (electrospray ioniza- 
tion) source (Thermo Finnigan) and HP 1100 liquid chro- 
matography system (Agilent Technologies) in negative ion 
detection mode. Nucleoside analysis was performed as de- 
scribed previously (31). In brief, about 4 (xg of isolated 
tRNA was digested to nucleosides with nuclease PI (Wako 
Pure Chemical Industries) and bacterial alkaline phos- 
phatase C75 (Takara bio), and then analyzed on the LCQ 
DUO ion-trap mass spectrometer with an ESI source and 
HP 1100 liquid chromatography system in positive-ion de- 
tection mode. ProMass (Novatia) was used to obtain the 
uncharged masses of whole tRNAs by deconvolution of the 
multiply charged mass spectra. 

RESULTS 

Isolation of 22 species of bovine mt tRNAs 

To date, 1 1 species of bovine mt tRNAs have been described 
in published work, along with the details of their post- 
transcriptional modifications: Ser(UCN) (33), Ser(AGY) 
(34,35), Phe (36), Arg (35), Gly (35), He (35), Met (8,37), 
Val (35), Trp (35), Glu (7) and Gin (7). In addition, se- 
quences of four species of bovine mt tRNAs for Thr, 
Leu(UUR), Leu(CUN) and Lys have recently been de- 
posited in the tRNA databases as personal communications 
(38,39). However, some tRNA sequences deposited in the 



tRNA databases contain unidentified modified bases and 
mis-annotation of modifications determined by outmoded 
methods. To accurately define the complete landscape of 
post-transcriptional modifications of mammalian mt tR- 
NAs, we analyzed all 22 species of mt tRNAs, which we pre- 
viously isolated from bovine liver by chaplet column chro- 
matography (29). 

Mass spectrometric analysis of post-transcriptional modifi- 
cations in bovine mt tRNAs 

We performed a series of mass spectrometric analy- 
ses of each individual mt tRNA to determine its post- 
transcriptional modifications. As an example, determina- 
tion of post-transcriptional modifications in mt tRNA Ala 
is depicted in Figure 1A. First, nucleoside analysis by 
LC/ESI-MS was carried out to determine the composition 
of modifications, revealing that 1-methyladenosine (ir^A), 
TV^-methylguanosine (m 2 G) and pseudouridine (^) were 
present in mt tRNA Ala (Figure IB). In parallel, we used cap- 
illary LC/nano ESI-MS to analyze RNA fragments gen- 
erated by digestion with RNase Ti (Figure 1C) or RNase 
A (Figure ID). Unmodified fragments were identified by 
comparing the observed m/z values with the calculated 
m/z values deduced from the mt tRNA Ala gene encoded 
in bovine mtDNA (GenBank accession number V00654). 
Nucleotide compositions of internal fragments were calcu- 
lated as 5 / -hydroxyl (5'OH) and 3 / -phosphate (3'p). The se- 
quence of each fragment was further analyzed by collision- 
induced dissociation (CID) (Figure IE). Because tRNA 
has 5 / -phosphate and 3 / -hydroxyl groups, the 5 / -terminal 
fragment pGAGGAUp (MW 2097.2) and the 3' terminal 
fragment oh-CAAUCCUUACCA- oh (MW 3697.5) were 
detected as species with unique molecular masses distinct 
from those of the internal fragments. Furthermore, RNA 
fragments containing modified nucleosides could be iden- 
tified from the deduced molecular mass calculated from 
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Figure 1. Mass spectrometric analysis of bovine mt tRNA Ala for assignment of post-transcriptional modifications. (A) Secondary structure of bovine mi- 
tochondrial tRNA Ala with post-transcriptional modifications determined in this study. The position numbers of the modifications are displayed according 
to the nucleotide numbering system from the tRNA compilation (40). Symbols for modified nucleosides are as follows: n^A, 1 -methyladenosine; m 2 G, 
7V 2 -methylguanosine and ^, pseudouridine. Watson-Crick base pairs are indicated by solid lines, whereas G-U pairs are indicated by asterisks. (B) Nucle- 
oside analysis of bovine mitochondrial tRNA Ala . Left, top panel: UV chromatogram at 254 nm of the four major nucleosides (C, U, G and A). Left, lower 
panels: extracted-ion chromatograms (XIC) for the protonated ion of m 1 A nucleoside (m/z 282, black line) with its base ion (m/z 150, gray line) (second 
panel), m 2 G nucleoside (m/z 298, black line) with its base ion (m/z 166, gray line) (third panel) and ^ nucleoside (m/z 245, black line) (bottom panel). The 
XIC for the base ion (20% upper offset) is overlaid on the XIC for the nucleoside ion. Right: mass spectra of m 1 A and m 2 G Cleavage positions for the 
base-related ions are indicated on the chemical structures. (C) RNA fragment analysis of RNase Ti digests of bovine mitochondrial tRNA Ala . Assigned 
fragments are indicated on the base peak chromatogram (BPC) in the first panel, 'p' stands for the terminal phosphate group. The XIC for the doubly 
charged negative ion of a modification-containing fragment (AUUUm 1 Am 2 Gp, m/z 982.6) is indicated in the second panel. Because m 2 G at position 10 is 
a partial modification, both AUUUm 1 Am 2 Gp and AUUUm 1 AGp were detected. (D) RNA fragment analysis of RNase A digests of bovine mitochondrial 
tRNA Ala . Assigned fragments are indicated on the BPC. The XIC for the doubly charged negative ion of the 5'-terminal fragment (pGAGGAUp, m/z 
1047.6) is indicated in the second panel. (E) A CID spectrum of a cyanoethylated RNA fragment to determine the location of a ^ site. The doubly charged 
negative ion of the RNA fragment (m/z 1617.7) shown in the inset was used as the precursor ion for CID. The product ions were assigned according 
to McLuckey et al. (41). The asterisks in the spectrum denote product ions containing ce 1 ^. (F) Whole mass analysis of intact bovine mt tRNA Ala . A 
series of multiply charged negative ions is shown in the mass spectrum. The charge values are indicated in parentheses. The observed mass obtained by 
deconvoluting the mass spectra is shown in the inset. 
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the combinations of unmodified and expected modified 
nucleotides. For example, in the RNase Ti digest of mt 
tRNA Ala , we found a species with m/z 982.627, correspond- 
ing to a doubly charged negative ion of dimethylated AUU- 
UAGp. The positions of modifications were determined by 
CID analysis, which indicated the modified sequence AU- 
UUA(+me)G(+me)p. According to the nucleoside analysis, 
this species was assigned as AUUUm 1 Am 2 Gp, with n^A 
and m 2 G at positions 9 and 10, respectively. Similarly, m 2 G 
was also found at position 26. 

To identify ^, a mass-silent modification, each tRNA 
was treated with acrylonitrile to cyanoethylate ^ (1- 
cyanoethyl ^; ce 1 ^). The derivatized tRNA was then di- 
gested by RNase Ti and subjected to mass spectrometry to 
detect the cyanoethylated fragments with molecular mass 
increased by 53 Da. Three cyanoethylated fragments were 
detected in the RNase Ti digest, and each fragment was 
further analyzed by CID. As shown in Figure IE, ^ at posi- 
tion 39 was identified from the assignment of product ions 
of CAUce^CAAUUGp. Accordingly, we identified three 
^s at positions 28, 39 and 57. Finally, to confirm the assign- 
ment of post-transcriptional modifications, the total molec- 
ular mass of mt tRNA Ala was measured by deconvoluting 
the multiply charged negative ions produced by ESI (Fig- 
ure IF). The observed mass (23178.5 Da) was fairly close 
to the calculated mass (23180.6 Da), which corresponds to 
the base composition of mt tRNA Ala (PU27C9A20G16) with 
three methyl groups. Other tRNAs were basically analyzed 
using the same procedure as for mt tRNA Ala (Supplemen- 
tary Information and Table SI). 

Post-transcriptional modifications in seven bovine mt tRNAs 

For 10 of the species of mt tRNAs described in the 
previous literature, our determinations of modifications 
yielded results consistent with published findings: Phe, Gly, 
Leu(CUN), Met, Arg, Ser(UCN), Ser(AGY), Val, Glu and 
Gin (Supplementary Figure SI). 

Several of the tRNAs we characterized have not been pre- 
viously reported or deposited in any databases. We deter- 
mined the primary structures and modifications for seven 
such mt tRNAs: Cys, Asp, His, Asn, Pro and Tyr (Figure 
2 and Supplementary Table SI) plus the aforementioned 
mt tRNA Ala (Figure 1). We identified queuosine (Q) at the 
wobble positions of mt tRNAs for Asp, His, Asn and Tyr. 
In cytosolic tRNA Asp and tRNA Tyr , Q was glycosylated to 
mannosyl-Q and galactosyl-Q, respectively; however, Q was 
not glycosylated in the mt tRNAs. The mt tRNAs for Ala 
and Pro correspond to family boxes; as expected, unmodi- 
fied U was present at the wobble position. In mt tRNA Hls , 
we identified a G at position -1 (G_i) (Supplementary Fig- 
ure S2). This base is likely to be added enzymatically af- 
ter transcription and 5 / -end processing. In mt tRNA Asp , 
m 2 G was present at position 6, the first such case to be re- 
ported in mammalian mitochondria. At position 37, which 
is 3 / -adjacent to the anticodon, A^-isopentenyladenosine 
(i 6 A) and 2-methylthio-7V^-isopentenyladenosine (ms 2 i 6 A) 
were present in mt tRNA Cys and mt tRNA Tyr , respectively. 
A^-threonylcarbamoyladenosine (t 6 A) and n^G were also 
present in the mt tRNAs for Asn and Pro, respectively. Dur- 
ing our analysis, we found two polymorphisms in the se- 



quences of mt tRNAs for Asp and Tyr: in mt tRNA Asp , a 
mixture of A and G were present at position 57; similarly, A 
and G were mixed at position 5 in mt tRNA Tyr (Figure 2). 
Both sites are encoded as A7356 (H-chain) and A5681 (L- 
chain) in the bovine mtDNA sequence (GenBank V00654) 
used as a reference. We speculate that these polymorphisms 
are the result of mtDNA heteroplasmy in the bovine liver 
we analyzed. In fact, both sites are encoded as guanosines in 
the mtDNA sequence of Bos taurus isolate FC3 (GenBank 
accession number DQ124389.1), suggesting that these sites 
are polymorphic in populations of healthy individuals. 



Full complement of modifications in five bovine mt tRNAs 

Although sequences and modifications of the remaining 
five mt tRNAs [He, Leu(UUR), Lys, Thr and Trp] have 
been published or deposited in the database (38,39), we dis- 
covered eight previously unreported modifications in these 
tRNAs (Figure 3): TV 2 , A^ 2 -dimethylguanosine (m 2 2 G) and 
m 1 A at positions 26 and 58 in mt tRNA Ile ; Tm 5 U at position 
34 in mt tRNAs for Leu(UUR) and Trp; Tm 5 s 2 U and t 6 A 
at positions 34 and 37 in mt tRNA Lys and 5-methylcytidine 
(m 5 C) at position 72 in mt tRNA Thr and position 48 in mt 
tRNA Trp . The wobble modifications of bovine mt tRNAs 
for Leu(UUR), Lys and Trp remained unidentified in tR- 
NAdb (38). In a previous study (6), we reported that im 5 U 
and Tm 5 s 2 U were present in bovine mt tRNA Leu ^ UUR ^ and 
mt tRNA Lys , respectively. To date, however, their full se- 
quences had not been reported. The m 5 C72 in mt tRNA Thr 
is the first reported instance of this modification in mam- 
malian mitochondria. 



DISCUSSION 

Here, we report the post-transcriptional modifications in 
seven species of bovine mt tRNAs not previously character- 
ized, and eight previously unidentified modified bases in five 
mt tRNAs whose sequences and modifications were deter- 
mined in earlier studies. In total, we identified 1 5 species of 
modified nucleosides at 1 18 positions in the complete set of 
bovine mt tRNAs (Supplementary Table S2), i.e. 7.48% of 
the bases in these tRNAs are modified. The sites and species 
of all modifications are summarized in a schematic clover- 
leaf structure (Figure 4). Notably, all modifications are base 
modifications; the absence of ^-O-methylation is a charac- 
teristic feature of mt tRNAs. 

Eight mt tRNAs responsible for family boxes had un- 
modified U at their wobble positions (Table 1), suggesting 
that the family boxes in mitochondria are decoded by sin- 
gle tRNAs via the four- way wobble rule. For NNR codons, 
all six mt tRNAs had wobble modifications (Table 1): f^C 
in mt tRNA Met , im 5 U in mt tRNAs for Leu(UUR) and 
Trp and Tm 5 s 2 U in mt tRNAs for Lys, Glu and Gin. We 
consider that these modifications are required for efficient 
recognition of the cognate codons, as well as to prevent mis- 
reading of near-cognate codons. Q was present at the wob- 
ble position in four mt tRNAs (Tyr, His, Asn and Asp) re- 
sponsible for NAY codons (Table 1). Q34 is known to re- 
strict the conformational flexibility of the anticodon loop 
by making hydrogen bonds between the side chain amine 
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Figure 2. Post-transcriptional modifications in six bovine mt tRNAs. Symbols for modified nucleosides are as follows: n^G, 1-methylguanosine; i 6 A, 
TV^-isopentenyladenosine; Q, queuosine; m 5 C, 5-methylcytidine; ms 2 i 6 A, 2-methylthio-7V^-isopentenyladenosine. The 'G/A' in tRNA Asp and tRNA Tyr 
indicates that both G and A were detected at this position, probably as a result of heteroplasmy in the bovine liver we used in this study 



of Q-base and 2' OH of U33 (42). Q34 enables tRNA to de- 
code NAU efficiently (43). Although the functional role of 
Q still remains obscure, it is associated with various phys- 
iological events in cytoplasmic tRNAs, including cell de- 
velopment and proliferation, neoplastic transformation and 
translational read-through or frameshift essential for retro- 
viral production (44^6). These facts suggest that the pres- 
ence of Q in the four mt tRNAs plays a modulatory role in 
deciphering NAY codons. 

Other mt tRNAs responsible for NNY codons [Phe, He, 
Cys, Ser(AGY)] had unmodified G at the wobble position 
(Table 1). The results clearly reveal a general principle to 
decipher the minimal decoding system by the base modifi- 
cations at the wobble positions. 

Hypermodifications are frequently found at position 37 
of tRNAs. These modifications play important roles in 
maintaining translational efficiency and integrity (47,48). 
Four species of base modifications were present at position 
37 in 13 mt tRNAs (Figure 4 and Supplementary Figure 
SI). Among six mt tRNAs responsible for UNN codons, 
five had i 6 A or ms 2 i 6 A at position 37. t 6 A37 was present in 
the mt tRNAs for He, Thr, Asn, Lys and Ser(AGY), whereas 
m 1 G37 was present in the mt tRNAs for Leu(CUN), Pro 
and Gin. Last year, we discovered cyclic t 6 A, which is 
formed by ATP-dependent dehydration of t 6 A catalyzed by 
TcdA, as a bonafide modification at position 37 of tRNAs 
from bacteria, yeast, plants and protists (49). However, we 



did not detect any cyclic form of t 6 A in bovine mt tRNAs, 
and there is no homolog of TcdA in mammalian genomes. 

Loss-of-function mutations in tRNA-modifying en- 
zymes can cause human diseases. Consistent with this idea, 
large-scale disease-associated exome analyses have identi- 
fied a number of genes that encode tRNA-modifying en- 
zymes (50,51). Identification of all genes involved in mt 
tRNA modifications will help us to identify genes and mu- 
tations associated with various diseases, especially those 
linked to mitochondrial dysfunction. Basically, the struc- 
ture and sequence of human mt tRNAs are similar to those 
of bovine mt tRNAs. In fact, 81% of total bases in 22 mt 
tRNAs are conserved between these two mammals. That is 
why human mt tRNA-modifying enzymes can be predicted 
based on information of bovine mt tRNAs. So far, nine hu- 
man genes have been confirmed to be responsible for base 
modifications in mammalian mt tRNAs, and others have 
been predicted based on studies of tRNA modifications in 
model organisms (Table 2). Several of these modifications, 
and the enzymes that may catalyze their formation, are de- 
scribed in the following paragraphs. 

On the basis of studies of bacterial and yeast mitochon- 
dria (63,78), GTPBP3 and MTOl are predicted to be en- 
zymes involved in Tm 5 U formation; however, to date, no 
direct evidence in support of this prediction has been pub- 
lished. The GTPBP3-MT01 complex has been proposed to 
recognize mt tRNAs for Leu(UUR), Trp, Lys, Glu and Glu. 
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Table 2. List of confirmed and predicted genes responsible for post-transcriptional modifications in mammalian mt tRNAs 



Confirmed gene(s) in 

Position a tRNA species Modification 13 human or mammals Predicted gene(s) in human 



6 


Asp 


m 2 G 




THUMPD3 or THUMPD2 










(52,53) 


9 


Ala, Asp, Glu, Phe, Gly, 


m'A 


TRMT10C and SDR5C1 






His, Lys, Leu(CUN),Asn, 




(54) 






Pro, Arg, Thr, Val, Trp 










Cys, He, Leu(UUR), Gin, 


m'G 


TRMT10C and SDR5C1 






y 




(54) 




10 


Ala, Phe, Gly, His, Lys, 


m 2 G 




TRMT11 and TRMT112 (55) 




Leu(UUR), Leu(CUN), 










Asn, Pro, Val, Trp, Tyr 








20 


LeuflJUR^ LeufCUN^ 


D 








Ser(UCN) 








OA 


Ala fl-1-n T mi^TITTP^ 

/Via, OIU, i^eu^U UrvJ 


m kj 




1 K1V1 1 1 




lie 


m z 2G 




TRMT1 (57,58) 


2 /a 


^er(UCJN) 








27 


A <5t"> T— Ti o Tip 

X11l5, Alt/, 




















A«?n Pro Val lVfpt 








28 


Ala, Cys, Glu, Lys, 




PUS1 (25) 






Leu'(CUN), Asn, 










Ser(UCN), Tyr 








29 


Ser(UCN) 






PUS1 


31 


Asp, Leu(CUN) 






RPUSD1, RPUSD2, 








RPUSD3 or RPUSD4 d 


32 


Cys, Val 


4> 




RPUSD1, RPUSD2, 








RPUSD3 or RPUSD4 d 




Ser(UCN), Thr 


m 3 C 




METTL2B (59,60) 


34 


Leu(UUR), Trp, Glu, Lys, 


Tm 5 U 




GTPBP3 and MTOl (61,62) 




Gin, 










Glu, Lys, Gin 


Tm 5 s 2 U e 


MTUP (63) and TW^ye 










(64) 






Met 


f 5 C 




Unidentified 




Asp, His, Asn, Tyr 


Q 




hQTRTl and hQTRTDl 










(65,66) 


37 


lie, Lys, Asn, Ser(AGY), 


t 6 A 




mDC and 07?/7 (OSGEPL1) 




Thr 






(67) 




Cys, Phe, Ser(UCN), Trp, 


i 6 A 


TRIT1 (68) 






Tyr 










Phe, Ser(UCN), Trp, Tyr 


ms 2 i 6 A f 


CDK5RAP1^ (69) 






Leu(CUN), Pro, Gin 


m'G 


(70) 




39 


Ala, Cys, Phe, Gly, His, 






PraJ (71) 




Leu(U UK), Lrm, Arg, lyr 








40 


Glu, Gin 






PtAS^ 


48 


Leu(UUR), Asn, Trp 


m 5 C 




(72,73) 


49 


Glu, Ser(AGY) 


m 5 C 




7VS£/7V2 (72,73) 


50 


Met 


*I> 




Unidentified 


55 


Glu, Gin, Ser(UCN), Tyr 


4> 




7W£2 (74) 


57 


Ala 






Unidentified 


58 


Cys, Glu, lie, Lys, 


m'A 


TRMT61B (75) 






Leu(UUR), Ser(UCN) 








67 


Thr 






P*AS7 d (76) 


72 


Thr 


m 5 C 




Unidentified 



a Numbering system for tRNA comes from the tRNAdb compilation (40). 

b Symbols for modifications originate from MODOMICS (http://modomics.genesilico.pl/) (39). 

c This position number is unique to mt tRNA Ser(UCN) (77). 

d These predictions were altered from our previous prediction (2). 

e MTUl and NFS1 are involved in 2-thiolation of Tm 5 s 2 U34. MTU1 (Mitochondrial tRNA-specific 2-thiouridylase 1) is known as TRMU which originated 
from bacterial trmll (tRNA methyltransferase U). However, trmU (renamed as mnmA) was found to be a mis-annotation because it is not a tRNA 
methyltransferase. 

f CDK5RAPl is required for 2-methylthiolation of ms 2 i 6 A37. 
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Figure 3. Revised information regarding post-transcriptional modifications of five bovine mt tRNAs. The updated modified bases are represented in bold 
type. Symbols for modified nucleoside are as follows: m 2 2G, N 2 , 7V 2 -dimethylguanosine; D, dihydrouridine; Tm 5 U, 5-taurinomethyluridine; Tm 5 s 2 U, 5- 
taurinomethyl-2-thiouridine; t 6 A, TV^-threonylcarbamoyladenosine; m 3 C, 3-methylcytidine. 



However, these five mt tRNAs share little sequence simi- 
larity and no consensus motif. Nonetheless, Tm 5 U forma- 
tion is sensitive to a single pathogenic point mutation asso- 
ciated with MEL AS or MERRF (12), implying strict sub- 
strate specificity of the enzyme. To resolve this issue, it will 
be necessary to perform in vitro reconstitution of im 5 U on 
mt tRNAs using recombinant proteins. 

We previously identified MTU1, a thiouridylase respon- 
sible for 2-thiolation of Tm 5 s 2 U in mt tRNAs (63). The sul- 
fur atom comes from Cys, a process mediated by a cysteine 
desulfurase, NFS1 (64). During bacterial 2-thiouridine for- 
mation, several sulfur mediators transfer persulfide sulfur 
from cysteine desulfurase to thiouridylase (79). To date, 
however, there is no information available regarding the in- 
volvement of such sulfur mediators in the formation of 2- 
thiouridine in mitochondria. 

In mammals, the Q-base is obtained either from the diet 
or through intestinal microflora (45), because no homolo- 
gous genes responsible for Q-base biogenesis de novo are 
encoded in mammalian genome. Mammals possess two 
tRNA-guanine transglycosylases (TGTases), QTRT1 and 
QTRTD 1 , which may transglycosylate the wobble base with 
the dietary Q-base. According to its subcellular localization, 
QTRTD1 appears to be the mitochondrial TGTase (65,66). 
Further studies will be necessary to confirm the identity of 
the mitochondrial TGTase, as well as the transport path- 
way by which dietary Q enters the mitochondria. Q is not 
present in mt tRNA Asp from Morris hepatoma cells (80). 



In yeast, t 6 A37 is synthesized by two enzymes, Sua5p and 
Qri7p (67). YRDC and OSGEPL1, the human homologs 
of yeast Sua5p and Qri7p, are presumably involved in t 6 A 
formation in mammalian mt tRNAs. 

TRIT1 is a tumor suppressor gene, mutations in which 
are associated with cancer progression. Human TRIT1 was 
identified as a tRNA isopentenyltransferase (IPTase) ho- 
mologous to bacterial MiaA and yeast Mod5; the latter is 
responsible for i 6 A37 formation in both cytoplasmic and mt 
tRNAs. Knockdown of TRIT1 in human cells reduced the 
abundance of i 6 A in mt tRNA Ser(UCN ) (68). In the present 
study, we determined that bovine mt tRNAs for Tyr and 
Cys have i 6 A37; in total, five mt tRNAs have been identi- 
fied as potential substrates for TRIT1 in mitochondria. 

The human protein CDK5RAP1 is a 2-methylthiolase 
that is homologous to bacterial MiaB. A CDK5RAP1 
variant localizes in mitochondria, and knockdown of 
CDK5RAP1 results in a reduction of ms 2 i 6 A levels (69). 
Thus, four mt tRNAs [Phe, Ser(UCN), Tyr and Trp] are 
likely to be modified by CDK5RAP1. 

TRMT5, a methyltransferase responsible for m 1 G37 for- 
mation, is predicted to be localized to the mitochondria 
(70,81), and recombinant TRMT5 can modify in vitro tran- 
scribed mt tRNA (70), suggesting that TRMT5 methylates 
three mt tRNAs: Leu(CUN), Pro and Gin. 

Base methylation at position 9 is frequently observed 
among bovine mt tRNAs. Indeed, 19 mt tRNAs have m 1 A9 
or m 1 G9 (Figure 4). m^9 stabilizes the canonical clover- 
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Figure 4. Summary of post-transcriptional modifications in bovine mt tR- 
NAs. Species and numbers of post-transcriptional modifications identified 
in 22 bovine mt tRNAs are mapped on the schematic secondary structure 
of tRNA. The modified positions are depicted by gray circles with a symbol 
indicating each modification. At each position, the number of tRNAs that 
bear the modification are shown in parenthesis. Positions 27a and 43a, in- 
dicated by light gray circles, are unique to mt tRNA Ser ( UCN ) . G_i is specific 
to mt tRNA His . 



leaf structure of mt tRNAs (82). In addition, m 1 A9 is in- 
dispensable for aminoacylation of nematode mt tRNAs 
that lack the T-arm. The methyltransferase for m^9 or 
m 1 G9 has been identified as a complex of TRM10C and 
SDR5C1, both of which are components of mt RNase P 
(54). The TRM10C-SDR5C1 complex has broad substrate 
specificity for 19 out of 22 mt tRNA species. Only three 
mt tRNAs are unmethylated at position 9. Of these, mt 
tRNA Met has C at this position, and the two mt tRNA Ser 
isoacceptors have non-canonical cl overleaf structures: Type 
I for mt tRNA Ser ( UCN ) and Type III for mt tRNA Ser < AGY ) (2). 

^ is the most abundant modification in bovine mt tR- 
NAs. Among the 22 mt tRNAs, we mapped ^ at 42 sites, 
35 of which reside in the anticodon arm. In humans, there 
are 13 species of ^ synthases, each of which has different 
substrate specificity and introduces ^ at multiple sites in 
various RNAs. PUS1 is a human synthase with broad 
substrate specificity that introduces ^ into cytoplasmic tR- 
NAs, mt tRNAs, U2 snRNA and SRA non-coding RNA 
(25,76,83-84), and it is responsible for ^27 and ^28 in mt 
tRNAs (25). We assume that PUS1 is also responsible for 
^I>27a and ^29 in mt tRNA Ser( ^ UCN ^ and ^67 in mt tRNA Thr 
(76). Because missense mutations in the PUS1 gene are re- 
sponsible for the development of MLASA (25), it will be 
necessary to determine the ^ sites targeted by PUS1 to gain 
a deeper understanding of the functional role of this modi- 
fication, as well as the molecular pathogenesis of MLASA. 



According to the studies of yeast and Escherichia coli tR- 
NAs, ^39 and ^40 are likely to be introduced by PUS3 (71), 
whereas ^55 is probably modified by TRUB2 (74). In yeast 
mitochondria, ^31 and ^32 are introduced by PUS6 and 
PUS9, respectively (85,86). Although we previously pre- 
dicted that either RPUSD2 or RPUSD4 was responsible 
for ^31 and ^32 (2), two other homologs, RPUSD1 and 
RPUSD3, should also be regarded as candidates for these 
modifications. Regarding ^50 and ^57, no candidate en- 
zymes can be predicted at present. Biochemical and genetic 
analyses of ^ synthase genes will be necessary to determine 
the enzyme responsible for each ^ site. 

In contrast to cytoplasmic tRNAs, most mammalian mt 
tRNAs do not have consensus sequences of D- and T- 
loops, which are conserved among tRNAs in general; con- 
sequently, they have lost the canonical D-loop/T-loop inter- 
action. This feature is characteristic of mt tRNAs in mam- 
mals (2). In the D-loop, dihydrouridines (D) at positions 
16 and 17, which are frequently observed in canonical tR- 
NAs from various sources, are not present in bovine mt tR- 
NAs. Instead, D20 was present in three species of bovine 
mt tRNAs. Based on a study on yeast tRNA, it is reason- 
able to predict that human tRNA-dihydrouridine synthase 
2 (DUS2) introduces D20 (56). Human DUS2 has been im- 
plicated in pulmonary carcinogenesis (87). In the T-loop, 
^55 and m 1 A58, both of which are typical modifications 
in canonical tRNAs, are present in four and six mt tRNAs, 
respectively. We previously identified human TRMT61B as 
the methyltransferase responsible for m 1 A58 in mt tRNAs 
for Leu(UUR), Ser(UCN) and Lys (75). In addition to 
these three tRNAs, we also found m 1 A58 in mt tRNAs for 
Cys, Glu and He. 5-methyluridine (m 5 U) at position 54 is 
one of the most common T-loop modifications in canoni- 
cal tRNAs. m 5 U54 is present in the human mt tRNAs for 
Leu(UUR) (19) and Ser(UCN) (88); however, we did not 
detect m 5 U54 in any bovine mt tRNAs. 

m 2 G6 is present in bovine mt tRNA Asp . Human mt 
tRNA Asp does not have m 2 G6 (89), because this position 
is replaced by A6. The methyltransferase responsible for 
m 2 G6 was identified as Trml4 in Methanocaldo coccus jan- 
naschii (52) and TTH1157 in Thermus thermophilus (53). A 
similarity search using the Trml4 sequence retrieved two 
human homologs, THUMPD2 and THUMPD3. Accord- 
ing to the WoLF P-sort (90) prediction of subcellular lo- 
calization, the former is likely localized in the cytoplasm 
and the latter in mitochondria. Although further investiga- 
tion is necessary, THUMPD3 is the most plausible candi- 
date for the enzyme that introduces m 2 G6 in mt tRNA Asp , 
whereas THUMPD2 probably acts on cytoplasmic tRNAs. 
Although both genes are encoded in the human genome, 
no m 2 G6 is present in human mt tRNA Asp , implying the 
presence of other substrates for this enzyme in human mi- 
tochondria. 

In this study, we identified G_i at the 5' terminus of 
bovine mt tRNA Hls . In general, G_i is an identity element 
for aminoacylation by histidyl-tRNA synthetase (91). This 
is the first reported instance of this base in mammalian 
mitochondria, although G_i is present in both yeast and 
starfish mt tRNA His (92,93). Judging from the fact that the 
5 / -adjacent nucleotide of the mt tRNA Hls gene in bovine 
mt DNA is T rather than G (e.g. GenBank V00654), it is 
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likely that G_i is added enzymatically after S^terminal pro- 
cessing by mt RNase P. THG1L is a mammalian guany- 
ly transferase that adds G_i at the 5' terminus of cyto- 
plasmic tRNA Hls . Because THG1L is predicted to have a 
mitochondria-targeting sequence at the N- terminus (94), 
THG1L is likely to be responsible for G_i addition of mt 
tRNA His . 

In summary, we have compiled the first complete picture 
of post-transcriptional modifications in mammalian mt tR- 
NAs. The results of this study enable a deeper understand- 
ing of the molecular mechanisms underpinning the minimal 
decoding system in mammalian mitochondria, and should 
help predict the human tRNA-modifying enzymes respon- 
sible for each modification in mt tRNAs. The list of tRNA- 
modifying enzymes serves as a practical landmark that en- 
courages us to identify all genes responsible for tRNA mod- 
ifications in mammalian mitochondria, and to further in- 
vestigate human diseases caused by tRNA modification dis- 
orders in mitochondria. An important goal for future ef- 
forts is the identification of all modifications in human mt 
tRNAs. 

SUPPLEMENTARY DATA 

Supplementary Data are available at NAR Online. 



ACKNOWLEDGEMENTS 

We are grateful to the members of the Suzuki laboratory, 
in particular, Dr. K. Miyauchi and Dr. A. Nagao, for many 
insightful discussions. We thank Prof. K. Watanabe (Tokyo 
University of Pharmacy and Life Sciences) for his abiding 
encouragement for our mitochondrial tRNA work. 

FUNDING 

Grants-in-Aid for Scientific Research on Priority Areas 
from the Ministry of Education, Science, Sports, and Cul- 
ture of Japan [to Takeo S. and Tsutomu S.]; JSPS Fellowship 
for Japanese Junior Scientists [to Takeo S.]. 
Conflict of interest statement. None declared. 



REFERENCES 

1. Wallace,D.C. (2005) A mitochondrial paradigm of metabolic and 
degenerative diseases, aging, and cancer: a dawn for evolutionary 
medicine. Annu. Rev. Genet., 39, 359-407. 

2. Suzuki,T., Nagao,A. and Suzuki,T. (2011) Human mtochondrial 
tRNAs: biogenesis, function, structural aspects, and diseases. Annu. 
Rev. Genet., 45, 299-329. 

3. Crick,F.H. (1966) Codon-anticodon pairing: the wobble hypothesis. 
J. Mol. Biol, 19, 548-555. 

4. Rogalski,M., Karcher,D. and Bock,R. (2008) Superwobbling 
facilitates translation with reduced tRNA sets. Nat. Struct. Mol. 
Biol, 15, 192-198. 

5. Yokoyama,S., Watanabe,T., Murao,K., Ishikura,H., Yamaizumi,Z., 
Nishimura,S. and Miyazawa,T. (1985) Molecular mechanism of 
codon recognition by tRNA species with modified uridine in the first 
position of the anticodon. Proc. Natl. Acad. Sci. U. S. A., 82, 
4905-4909. 

6. Suzuki,T., Suzuki,T., Wada,T., Saigo,K. and Watanabe,K. (2002) 
Taurine as a constituent of mitochondrial tRNAs: new insights into 
the functions of taurine and human mitochondrial diseases. EMBO 
J., 21, 6581-6589. 



7. Nagao,A., Suzuki,T., Katoh,T., Sakaguchi,Y. and Suzuki,T. (2009) 
Biogenesis of glutaminyl-mt tRNAGln in human mitochondria. Proc. 
Natl. Acad. Sci. U. S. A., 106, 16209-16214. 

8. Moriya,!, Yokogawa,T., Wakita,K., Ueda,T., Nishikawa,K., 
Crain,P.E, Hashizume,T., Pomerantz,S.C, McCloskeyJ.A., 
Kawai,G. et al. (1994) A novel modified nucleoside found at the first 
position of the anticodon of methionine tRNA from bovine liver 
mitochondria. Biochemistry, 33, 2234-2239. 

9. Takemoto,C, Spremulli,L.L., Benkowski,L.A., Ueda,T., 
Yokogawa,T. and Watanabe,K. (2009) Unconventional decoding of 
the AUA codon as methionine by mitochondrial tRNAMet with the 
anticodon f5CAU as revealed with a mitochondrial in vitro 
translation system. Nucleic Acids Res., 37, 1616-1627. 

10. Bilbille,Y, Gustilo,E.M., Harris,K.A., Jones,C.N, Lusic,H., 
Kaiser,R.I, Delaney,M.O., Spremulli,L.L., Deiters,A. and Agris,P.F. 
(201 1) The human mitochondrial tRNAMet: structure/function 
relationship of a unique modification in the decoding of 
unconventional codons. J. Mol. Biol. , 406, 251-214. 

11. Cantara,W.A., Murphy,F.V., Demirci,H. and Agris,P.F. (2013) 
Expanded use of sense codons is regulated by modified cytidines in 
tRNA. Proc. Natl. Acad. Sci. U.S.A., 110, 10964-10969. 

12. Suzuki,T, Nagao,A. and Suzuki,T. (2011) Human mitochondrial 
diseases caused by lack of taurine modification in mitochondrial 
tRNAs. Wiley Interdiscip Rev. RNA, 2, 376-386. 

13. Ruiz-Pesini,E., Lott,M.T, Procaccio,V., PooleJ.C, Brandon,M.C, 
Mishmar,D., Yi,C, Kreuziger,!, Baldi,P. and Wallace,D.C. (2007) 
An enhanced MITOMAP with a global mtDNA mutational 
phylogeny. Nucleic Acids Res., 35, D823-D828. 

14. Levinger,L., Morl,M. and Florentz,C. (2004) Mitochondrial tRNA 3' 
end metabolism and human disease. Nucleic Acids Res., 32, 
5430-5441. 

15. Belostotsky,R., Frishberg,Y. and Entelis,N. (2012) Human 
mitochondrial tRNA quality control in health and disease: a 
channelling mechanism? RNA Biol, 9, 33-39. 

16. Goto,Y, Nonaka,I. and Horai,S. (1990) A mutation in the 
tRNA(Leu)(UUR) gene associated with the MELAS subgroup of 
mitochondrial encephalomyopathies. Nature, 348, 651-653. 

17. Goto,Y, Nonaka,I. and Horai,S. (1991) A new mtDNA mutation 
associated with mitochondrial myopathy, encephalopathy, lactic 
acidosis and stroke-like episodes (MELAS). Biochim. Biophys. Acta, 
1097, 238-240. 

18. ShoffnerJ.M., Lott,M.T, Lezza,A.M., Seibel,R, Ballinger,S.W. and 
Wallace,D.C. (1990) Myoclonic epilepsy and ragged-red fiber disease 
(MERRF) is associated with a mitochondrial DNA tRNA(Lys) 
mutation. Cell, 61, 931-937. 

19. Yasukawa,T, Suzuki,T, Suzuki,T, Ueda,T, Ohta,S. and 
Watanabe,K. (2000) Modification defect at anticodon wobble 
nucleotide of mitochondrial tRNAs Leu (UUR) with pathogenic 
mutations of mitochondrial myopathy, encephalopathy, lactic 
acidosis, and stroke-like episodes. J. Biol. Chem., 275, 4251-4257. 

20. Yasukawa,T, Suzuki,T, Ishii,N, Ueda,T, Ohta,S. and Watanabe,K. 
(2000) Defect in modification at the anticodon wobble nucleotide of 
mitochondrial tRNA Lys with the MERRF encephalomyopathy 
pathogenic mutation. FEBS Lett, 467, 175-178. 

21. Kirino,Y, Goto,Y, Campos,Y, Arenas,! and Suzuki,T. (2005) 
Specific correlation between the wobble modification deficiency in 
mutant tRNAs and the clinical features of a human mitochondrial 
disease. Proc. Natl. Acad. Sci. U.S.A., 102, 7127-7132. 

22. Yasukawa,T, Kirino,Y, Ishii,N, HoltJ.J., Jacobs,H.T, 
Makifuchi,T, Fukuhara,N, Ohta,S., Suzuki,T. and Watanabe,K. 
(2005) Wobble modification deficiency in mutant tRNAs in patients 
with mitochondrial diseases. FEBS Lett. , 579, 2948-2952. 

23. Kirino,Y, Yasukawa,T, Ohta,S., Akira,S., Ishihara,K., Watanabe,K. 
and Suzuki,T. (2004) Codon-specific translational defect caused by a 
wobble modification deficiency in mutant tRNA from a human 
mitochondrial disease. Proc. Natl. Acad. Sci. U. S. A., 101, 
15070-15075. 

24. Yasukawa,T, Suzuki,T, Ishii,N, Ohta,S. and Watanabe,K. (2001) 
Wobble modification defect in tRNA disturbs codon-anticodon 
interaction in a mitochondrial disease. EMBO J., 20, 4794-4802. 

25. PattonJ.R., Bykhovskaya,Y, Mengesha,E., Bertolotto,C. and 
Fischel-Ghodsian,N. (2005) Mitochondrial myopathy and 
sideroblastic anemia (MLASA): missense mutation in the 



7356 Nucleic Acids Research, 2014, Vol 42, No. 11 



pseudouridine synthase 1 (PUS1) gene is associated with the loss of 
tRNA pseudouridylation. J. Biol. Chem., 280, 19823-19828. 

26. Zeharia,A., Shaag,A., Pappo,0., Mager-Heckel,A.M., Saada,A., 
Beinat,M., Karicheva,0., Mandel,H., Ofek,N., Segel,R. et al. (2009) 
Acute infantile liver failure due to mutations in the TRMU gene. Am. 
J. Hum. Genet., 85, 401-407. 

27. Ghezzi,D., Baruffini,E., Haack,T.B., Invernizzi,R, Melchionda,L., 
Dallabona,C, Strom,T.M., Parini,R., Burlina,A.B., Meitinger,T. 

et al. (2012) Mutations of the mitochondrial-tRNA modifier MTOl 
cause hypertrophic cardiomyopathy and lactic acidosis. Am. J. Hum. 
Genet., 90, 1079-1087. 

28. Baruffini,E., Dallabona,C, Invernizzi,R, YarhamJ.W, 
Melchionda,L., Blakely,E.L., Lamantea,E., Donnini,C, Santra,S., 
Vijayaraghavan,S. et al. (2013) MTOl mutations are associated with 
hypertrophic cardiomyopathy and lactic acidosis and cause 
respiratory chain deficiency in humans and yeast. Hum. Mutat. , 34, 
1501-1509. 

29. Suzuki,T. and Suzuki,T. (2007) Chaplet column chromatography: 
isolation of a large set of individual RNAs in a single step. Methods 
Enzymol, 425, 231-239. 

30. Mengel-Jorgensen,J. and Kirpekar,F. (2002) Detection of 
pseudouridine and other modifications in tRNA by cyanoethylation 
and MALDI mass spectrometry. Nucleic Acids Res., 30, el 35. 

31. Suzuki,T., Ikeuchi,Y, Noma,A., Suzuki,T. and Sakaguchi,Y. (2007) 
Mass spectrometric identification and characterization of 
RNA-modifying enzymes. Methods Enzymol, 425, 21 1-229. 

32. Soma,A., Ikeuchi,Y, Kanemasa,S., Kobayashi,K., Ogasawara,N., 
Ote,T., Kato,I, Watanabe,K., Sekine,Y. and Suzuki,T. (2003) An 
RNA-modifying enzyme that governs both the codon and amino acid 
specificities of isoleucine tRNA. Mol. Cell, 12, 689-698. 

33. Yokogawa,T., Watanabe,Y, Kumazawa,Y, Ueda,T., Hirao,L, 
Miura,K. and Watanabe,K. (1991) A novel cloverleaf structure found 
in mammalian mitochondrial tRNA(Ser) (UCN). Nucleic Acids Res., 
19,6101-6105. 

34. Arcari,P. and Brownlee,G.G. (1980) The nucleotide sequence of a 
small (3S) seryl-tRNA (anticodon GCU) from beef heart 
mitochondria. Nucleic Acids Res., 8, 5207-5212. 

35. Roe,B.A., WongJ.F.H., Chen,E.Y. and Armstrong,P.W. (1981) 
Sequence analysis of mammalian mitochondrial tRNAs in 
recombinant DNA. In: Walton, A.G (ed). Recombinant DNA: 
Proceedings of the Third Cleveland Symposium on Macromolecules, 
Cleveland, Ohio, 22-26 June 1981. pp. 167-176. 

36. Wakita,K., Watanabe,Y, Yokogawa,T, Kumazawa,Y, Nakamura,S., 
Ueda,T, Watanabe,K. and Nishikawa,K. (1994) Higher-order 
structure of bovine mitochondrial tRNA(Phe) lacking the 'conserved' 
GG and T psi CG sequences as inferred by enzymatic and chemical 
probing. Nucleic Acids Res., 22, 347-353. 

37. Dubin,D.T. and HsuChen,C.C. (1984) Sequence and structure of a 
methionine transfer RNA from mosquito mitochondria. Nucleic 
Acids Res., 12, 4185-4189. 

38. Juhling,E, Morl,M., Hartmann,R.K., Sprinzl,M., Stadler,P.F. and 
PutzJ. (2009) tRNAdb 2009: compilation of tRNA sequences and 
tRNA genes. Nucleic Acids Res., 37, D159-D162. 

39. Machnicka,M.A., Milanowska,K., Osman Oglou,0., Purta,E., 
Kurkowska,M., 01chowik,A., Januszewski,W., Kalinowski,S., 
Dunin-Horkawicz,S., Rother,K.M. et al. (2013) MODOMICS: a 
database of RNA modification pathways-2013 update. Nucleic Acids 
Res., 41, D262-D267. 

40. Sprinzl,M. and Vassilenko,K.S. (2005) Compilation of tRNA 
sequences and sequences of tRNA genes. Nucleic Acids Res., 33, 
D139-D140. 

41. Mcluckey,S.A., Vanberkel,G.J. and Glish,G.L. (1992) Tandem 
mass-spectrometry of small, multiply charged oligonucleotides. J. 
Am. Soc. Mass Spectromet. , 3, 60-70. 

42. Yokoyama,S., Miyazawa,T., Iitaka,Y, Yamaizumi,Z., Kasai,H. and 
Nishimura,S. (1979) Three-dimensional structure of hyper-modified 
nucleoside Q located in the wobbling position of tRNA. Nature, 282, 
107-109. 

43. Meier,E, Suter,B., Grosjean,H., Keith,G. and Kubli,E. (1985) 
Queuosine modification of the wobble base in tRNAHis influences 'in 
vivo' decoding properties. EMBO J., 4, 823-827. 

44. Iwata-Reuyl,D. and de Crecy-Lagard,V. (2009) Enzymatic formation 
of the 7-deazaguanosine hypermodified nucleosides oftRNA. In: 
Grosjean,H. (ed.), DNA and RNA Modification Enzymes: Structure, 



Mechanism, Function and Evolution. Landes Bioscience, Austin, 
TX, pp. 377-391. 

45. Vinayak,M. and Pathak,C. (2010) Queuosine modification oftRNA: 
its divergent role in cellular machinery. Biosci. Rep., 30, 135-148. 

46. Bjork,G. (1995) Biosynthesis and function of modified nucleosides. 
In: Soll,D. and Rajbhandary,U.L. (eds),t RNA: Structure, 
Biosynthesis and Function. American Society for 
Microbiology, Washington, DC, pp. 165-205. 

47. Suzuki,T. (2005) Biosynthesis and Function oftRNA Wobble 
Modifications. In: Grosjean,H. (ed.), Fine-Tuning of RNA 
Functions by Modification and Editing. Springer- Verlag. Berlin, 
Heidelberg, Vol. 12, pp. 23-69. 

48. El Yacoubi,B., Bailly,M. and de Crecy-Lagard,V. (2012) Biosynthesis 
and function of posttranscriptional modifications of transfer RNAs. 
Annu. Rev. Genet., 46, 69-95. 

49. Miyauchi,K., Kimura,S. and Suzuki,T. (2013) A cyclic form of 
N6-threonylcarbamoyladenosine as a widely distributed tRNA 
hypermodification. Nat. Chem. Biol, 9, 105-111. 

50. Bamshad,M.I, Ng,S.B., Bigham,A.W, Tabor,H.K., Emond,M.I, 
Nickerson,D.A. and Shendure,! (2011) Exome sequencing as a tool 
for Mendelian disease gene discovery. Nat. Rev. Genet., 12, 745-755. 

51. Torres,A.G, Batlle,E. and Ribas de Pouplana,L. (2014) Role of 
tRNA modifications in human diseases. Trends Mol. 
Med.,doi:10.1016/j.molmed.2014.01.008 

52. Menezes,S., Gaston,K.W, Krivos,K.L., Apolinario,E.E., Reich,N.O., 
Sowers, K.R., Limbach,P.A. and PeronaJ.J. (2011) Formation of 
m2G6 in Methanocaldococcus jannaschii tRNA catalyzed by the 
novel methyltransferase Trml4. Nucleic Acids Res., 39, 7641-7655. 

53. Roovers,M., Oudjama,Y, Fislage,M., BujnickiJ.M., Versees,W. and 
Droogmans,L. (2012) The open reading frame TTC1157 of Thermus 
thermophilus HB27 encodes the methyltransferase forming 
N(2)-methylguanosine at position 6 in tRNA. RNA, 18, 815-824. 

54. Vilardo,E., Nachbagauer,C, Buzet,A., Taschner,A., Holzmann,! 
and Rossmanith,W. (2012) A subcomplex of human mitochondrial 
RNase P is a bifunctional methyltransferase-extensive moonlighting 
in mitochondrial tRNA biogenesis. Nucleic Acids Res., 40, 
11583-11593. 

55. Purushothaman,S.K., BujnickiJ.M., Grosjean,H. and Lapeyre,B. 
(2005) Trml lp and Trml 12p are both required for the formation of 

2- methylguanosine at position 10 in yeast tRNA. Mol. Cell Biol, 25, 
4359-4370. 

56. Xing,F, Hiley,S.L., Hughes,T.R. and Phizicky,E.M. (2004) The 
specificities of four yeast dihydrouridine synthases for cytoplasmic 
tRNAs. J. Biol. Chem., 279, 17850-17860. 

57. Liu,J. and Straby,K.B. (2000) The human 
tRNA(m(2)(2)G(26))dimethyltransferase: functional expression and 
characterization of a cloned hTRMl gene. Nucleic Acids Res., 28, 
3445-3451. 

58. Ellis,S.R., Morales,M.J., LiJ.M., Hopper,A.K. and Martin,N.C. 
(1986) Isolation and characterization of the TRM1 locus, a gene 
essential for the N2,N2-dimethylguanosine modification of both 
mitochondrial and cytoplasmic tRNA in Saccharomyces cerevisiae. J. 
Biol. Chem., 261, 9703-9709. 

59. Noma,A., Yi,S., Katoh,T, Takai,Y, Suzuki,T. and Suzuki,T. (2011) 
Actin-binding protein ABP140 is a methyltransferase for 

3- methylcytidine at position 32 of tRNAs in Saccharomyces 
cerevisiae. RNA, 17, 1111-1119. 

60. D'Silva,S., Haider,S.J. and Phizicky,E.M. (2011) A domain of the 
actin binding protein Abpl40 is the yeast methyltransferase 
responsible for 3-methylcytidine modification in the tRNA 
anti-codon loop. RNA, 17, 1 100-1 1 10. 

61. Li,X. and Guan,M.X. (2002) A human mitochondrial GTP binding 
protein related to tRNA modification may modulate phenotypic 
expression of the deafness-associated mitochondrial 12S rRNA 
mutation. Mol. Cell Biol, 22, 7701-7711. 

62. Li,X., Li,R., Lin,X. and Guan,M.X. (2002) Isolation and 
characterization of the putative nuclear modifier gene MTOl 
involved in the pathogenesis of deafness-associated mitochondrial 12 
S rRNA A1555G mutation. J. Biol. Chem., 211, 27256-27264. 

63. Umeda,N., Suzuki,T, Yukawa,M., Ohya,Y, Shindo,H., Watanabe,K. 
and Suzuki,T. (2005) Mitochondria-specific RNA-modifying 
enzymes responsible for the biosynthesis of the wobble base in 
mitochondrial tRNAs. Implications for the molecular pathogenesis of 
human mitochondrial diseases. J. Biol. Chem., 280, 1613-1624. 



Nucleic Acids Research, 2014, Vol 42, No. 11 7357 



64. Nakai,Y, Umeda,N., Suzuki,T., Nakai,M., Hayashi,H., Watanabe,K. 
and Kagamiyama,H. (2004) Yeast Nfslp is involved in 
thio-modification of both mitochondrial and cytoplasmic tRNAs. J. 
Biol. Chem., 279, 12363-12368. 

65. Boland,C, Hayes,R, Santa-MariaJ., Nishimura,S. and Kelly, V. P. 
(2009) Queuosine formation in eukaryotic tRNA occurs via a 
mitochondria-localized heteromeric transglycosylase. J. Biol. Chem., 
284, 18218-18227. 

66. Chen,Y.C, Kelly, V.P., Stachura,S.V. and Garcia,G.A. (2010) 
Characterization of the human tRNA-guanine transglycosylase: 
confirmation of the heterodimeric subunit structure. RNA, 16, 
958-968. 

67. Wan,L.C, Mao,D.Y, Neculai,D., Strecker J., Chiovitti,D., 
KurinovJ., Poda,G., Thevakumaran,N., Yuan,R, Szilard,R.K. et al. 
(2013) Reconstitution and characterization of eukaryotic 
N6-threonylcarbamoylation of tRNA using a minimal enzyme 
system. Nucleic Acids Res., 41, 6332-6346. 

68. Lamichhane,T.N, Mattijssen,S. and Maraia,R.J. (2013) Human cells 
have a limited set of tRNA anticodon loop substrates of the tRNA 
isopentenyltransferase TRIT1 tumor suppressor. Mol. Cell Biol., 33, 
4900-4908. 

69. Reiter,V., Matschkal,D.M., Wagner,M., Globisch,D., 
Kneuttinger,A.C, Muller,M. and Carell,T. (2012) The CDK5 
repressor CDK5RAP1 is a methylthiotransferase acting on nuclear 
and mitochondrial RNA. Nucleic Aacids Res., 40, 6235-6240. 

70. Brule,H., Elliott,M., Redlak,M., Zehner,Z.E. and Holmes,W.M. 
(2004) Isolation and characterization of the human tRNA-(NlG37) 
methyltransferase (TRM5) and comparison to the Escherichia coli 
TrmD protein. Biochemistry, 43, 9243-9255. 

71. Lecointe,E, Simos,G, Sauer,A., Hurt,E.C, Motorin,Y. and 
Grosjean,H. (1998) Characterization of yeast protein Degl as 
pseudouridine synthase (Pus3) catalyzing the formation of psi 38 and 
psi 39 in tRNA anticodon loop. J. Biol. Chem., 273, 1316-1323. 

72. Brzezicha,B., Schmidt,M., MakalowskaJ., Jarmolowski,A., 
Pienkowska J. and Szweykowska-Kulinska,Z. (2006) Identification of 
human tRNA:m5C methyltransferase catalysing intron-dependent 
m5C formation in the first position of the anticodon of the pre-tRNA 
Leu (CAA). Nucleic Acids Res., 34, 6034-6043. 

73. Tuorto,F., Liebers,R., Musch,T., Schaefer,M., Hofmann,S., 
Kellner,S., Frye,M., Helm,M., Stoecklin,G. and Lyko,F. (2012) RNA 
cytosine methylation by Dnmt2 and NSun2 promotes tRNA stability 
and protein synthesis. Nat. Struct. Mol. Biol, 19, 900-905. 

74. Becker,H.E, Motorin,Y, Planta,R.J. and Grosjean,H. (1997) The 
yeast gene YNL292w encodes a pseudouridine synthase (Pus4) 
catalyzing the formation of psi55 in both mitochondrial and 
cytoplasmic tRNAs. Nucleic Acids Res., IS, 4493-4499. 

75. Chujo,T. and Suzuki,T. (2012) Trmt61B is a methyltransferase 
responsible for 1 -methyladenosine at position 58 of human 
mitochondrial tRNAs. RNA, 18, 2269-2276. 

76. Behm-Ansmant,L, Massenet, S., Immel,E, PattonJ.R., Motorin,Y. 
and Branlant,C. (2006) A previously unidentified activity of yeast and 
mouse RNA:pseudouridine synthases 1 (Puslp) on tRNAs. RNA, 12, 
1583-1593. 

77. Watanabe,Y, Kawai,G, Yokogawa,T, Hayashi,N, Kumazawa,Y, 
Ueda,T, Nishikawa,K., HiraoJ., Miura,K. and Watanabe,K. (1994) 
Higher-order structure of bovine mitochondrial tRNA(SerUGA): 
chemical modification and computer modeling. Nucleic Acids Res., 
22, 5378-5384. 

78. Armengod,M.E., MoukadiriJ., Prado,S., Ruiz-Partida,R., 
Benitez-Paez,A., Villarroya,M., Lomas,R., Garzon,M.I, 
Martinez-Zamora,A., Meseguer,S. et al. (2012) Enzymology of 



tRNA modification in the bacterial MnmEG pathway. Biochimie. , 94, 
1510-1520. 

79. Ikeuchi,Y, Shigi,N, Kato J., Nishimura,A. and Suzuki,T. (2006) 
Mechanistic insights into sulfur relay by multiple sulfur mediators 
involved in thiouridine biosynthesis at tRNA wobble positions. Mol. 
Cell, 21, 97-108. 

80. Randerath,E., Agrawal,H.P. and Randerath,K. (1984) Specific lack of 
the hypermodified nucleoside, queuosine, in hepatoma mitochondrial 
aspartate transfer RNA and its possible biological significance. 
Cancer Res., 44, 1167-1171. 

81. Lee,C, Kramer,G, Graham,D.E. and Appling,D.R. (2007) Yeast 
mitochondrial initiator tRNA is methylated at guanosine 37 by the 
Trm5-encoded tRNA (guanine-Nl-)-methyltransferase. J. Biol. 
Chem., 282, 27744-27753. 

82. Helm,M. (2006) Post-transcriptional nucleotide modification and 
alternative folding of RNA. Nucleic Acids Res., 34, 721-733. 

83. Zhao,X., PattonJ.R., Davis,S.L., Florence,B., Ames,S.J. and 
Spanjaard,R.A. (2004) Regulation of nuclear receptor activity by a 
pseudouridine synthase through posttranscriptional modification of 
steroid receptor RNA activator. Mol. Cell, 15, 549-558. 

84. Zhao,X., PattonJ.R., Ghosh,S.K., Fischel-Ghodsian,N, Shen,L. and 
Spanjaard,R.A. (2007) Pus3p- and Puslp-dependent 
pseudouridylation of steroid receptor RNA activator controls a 
functional switch that regulates nuclear receptor signaling. Mol. 
Endocrinol, 21, 686-699. 

85. AnsmantJ., Motorin,Y, Massenet,S., Grosjean,H. and Branlant,C. 
(2001) Identification and characterization of the tRNA:Psi 

31 -synthase (Pus6p) of Saccharomyces cerevisiae. J. Biol. Chem., 276, 
34934-34940. 

86. Behm- AnsmantJ., Grosjean,H., Massenet,S., Motorin,Y. and 
Branlant,C. (2004) Pseudouridylation at position 32 of mitochondrial 
and cytoplasmic tRNAs requires two distinct enzymes in 
Saccharomyces cerevisiae. J. Biol. Chem., 279, 52998-53006. 

87. Kato,T, Daigo,Y, Hayama,S., Ishikawa,N, Yamabuki,T, Ito,T, 
Miyamoto, M., Kondo,S. and Nakamura,Y. (2005) A novel human 
tRNA-dihydrouridine synthase involved in pulmonary 
carcinogenesis. Cancer Res., 65, 5638-5646. 

88. Toompuu,M., Yasukawa,T, Suzuki,T, Hakkinen,T, SpelbrinkJ.N, 
Watanabe,K. and Jacobs,H.T. (2002) The 7472insC mitochondrial 
DNA mutation impairs the synthesis and extent of aminoacylation of 
tRNASer(UCN) but not its structure or rate of turnover. J. Biol. 
Chem., 277, 22240-22250. 

89. Messmer,M., PutzJ., Suzuki,T, Suzuki,T, Sauter,C, Sissler,M. and 
Catherine,F. (2009) Tertiary network in mammalian mitochondrial 
tRNAAsp revealed by solution probing and phylogeny. Nucleic Acids 
Res., 37, 6881-6895. 

90. Horton,R, Park,K.I, Obayashi,T, Fujita,N, Harada,H., 
Adams-Collier,C.J. and Nakai,K. (2007) WoLF PSORT: protein 
localization predictor. Nucleic Acids Res., 35, W585-W587. 

91. Giege,R., Sissler,M. and Florentz,C. (1998) Universal rules and 
idiosyncratic features in tRNA identity. Nucleic Acids Res. , 26, 
5017-5035. 

92. Sibler,A.R, Martin,R.P. and Dirheimer,G. (1979) The nucleotide 
sequence of yeast mitochondrial histidine-tRNA. FEBS Lett. , 107, 
182-186. 

93. Tomita,K., Ueda,T. and Watanabe,K. (1999) The presence of 
pseudouridine in the anticodon alters the genetic code: a possible 
mechanism for assignment of the AAA lysine codon as asparagine in 
echinoderm mitochondria. Nucleic Acids Res., 27, 1683-1689. 

94. JackmanJ.E., Gott J.M. and Gray,M.W. (2012) Doing it in reverse: 
3'-to-5' polymerization by the Thgl superfamily. RNA, 18, 886-899. 



